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Quark spectral properties above T c from Dyson-Schwinger equations 
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We report on an analysis of the quark spectral representation at finite temperatures based on the 
quark propagator determined from its Dyson-Schwinger equation in Landau gauge. In Euclidean 
space we achieve nice agreement with recent results from quenched lattice QCD. We find different 
analytical properties of the quark propagator below and above the deconfinement transition. Using 
a variety of ansatze for the spectral function we then analyze the possible quasiparticle spectrum, in 
particular its quark mass and momentum dependence in the high temperature phase. This analysis 
is completed by an application of the Maximum Entropy Method, in principle allowing for any 
positive semi-definite spectral function. Our results motivate a more direct determination of the 
spectral function in the framework of Dyson-Schwinger equations. 

PACS numbers: 12.38.Mh,14.65.Bt,25.75.Nq 



I. INTRODUCTION 

The wealth of data produced by the Relativistic Heavy- 
Ion Collider (RHIC) has lead to new insights about the 
nature of the quark-gluon plasma (QGP), in particular in 
the strongly coupled regime just above the chiral phase 
transition (see Refs. [lHIj for recent reviews). Despite the 
strong coupling, phenomena like the constituent quark 
number scaling of elliptic flow [j| still suggest quarks 
as quasiparticle excitations. This should also be con- 
nected with the success of quasiparticle models in de- 
scribing thermodynamic properties in this regime of the 
QCD phase diagram [1,0] and also serves as a foundation 
for transport approaches like the one discussed in Ref . Q . 
In another context, dilepton production in a heavy ion 
collision has been related to the spectral properties of 
the thermalizcd quasiparticles and specifically to the dis- 
persion relation of quarks [9MTTI . We therefore conclude 
that a detailed understanding of the quasiparticle spec- 
trum, in particular close to the chiral phase transition, is 
desirable. 

Our quantity of interest is the quark spectral function, 
which encodes all information about the two-point cor- 
relation function, in particular also its poles and conse- 
quently the quark's dispersion relation and decay width. 
In the realm of weak coupling reliable results are obtained 
using the hard-thermal loop (HTL) expansion 
One finds the by now well-known pattern of two quasipar- 
ticle excitations and a continuum contribution stemming 
from a branch cut in the quark propagator due to Lan- 
dau damping, i.e. the absorption of a space-like quark by 
a hard gluon or hard antiquark. The two quasiparticles 
correspond to an ordinary quark with a positive ratio of 
chirality to helicity and to a collective mode dubbed 'plas- 
mino' with a negative chirality to helicity ratio. Both 
have thermal masses of order gT and decay widths of 
order g 2 T, where g is the coupling constant and T the 



temperature. 

Beyond systematic weak-coupling expansions, there is 
no straightforward approach for the determination of the 
spectral function or its characteristics. Some insight 
has been gained using model calculations, finding sim- 
ilar structures as in the weak coupling limit besides a 
possible additional third quasiparticle excitation which 
is only present at small momenta [l5|-[l3]- The width of 
the quasiparticle peaks, which is of order g 2 T in HTL 
approximation, might however become of the order of 
the thermal mass for couplings at the scale of the chi- 
ral phase transition [Lsl ]. This questions the existence of 
well-defined quasiparticles at least at small momenta. 

Ab initio calculations of correlators within lattice QCD 
are performed in Euclidean space and an analytic contin- 
uation that is necessary to determine the spectral func- 
tion is strictly speaking not possible. There are, however, 
approaches that aim to extract attributes of spectral 
functions from the numerical data. One of these is the 
Maximum Entropy Method (MEM) QjJEO], which has 
been applied in similar contexts, e.g. for extracting the 
spectral functions of mesons [111, [22[ • For cold and dense 
matter MEM has also been applied within the framework 
of Dyson-Schwinger equations (DSEs) in Ref. [23[. In 
some sense MEM corresponds to fitting the "most likely" 
spectral function to the data but without restricting the 
form of the spectral function. In contrast to this another 
method has been established in Refs. [HI HH- There 
the authors assumed a certain shape for the quark spec- 
tral function including a few fitting parameters which 
were determined from the numerical data. In such an 
approach one relies on physical guidance as e.g. given by 
the HTL approximation for the construction of a suitable 
ansatz for the spectral function. 

A complementary approach to the correlators of QCD 
are functional methods such as Dyson-Schwinger equa- 
tions and the functional rcnormalisation group. In re- 
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cent years, these methods have been successfully applied 
to problems such as the characterization of the chiral 
and deconfinement transition [26l - l32j |. or the properties 
of gluons in the high temperature phase j33rl35l ] . In this 
work we do the first steps towards the determination of 
spectral functions in the framework of Dyson-Schwinger 
equations beyond simple rainbow approximations. Par- 
ticularly due to our input for the gluon propagator the 
DSE calculations are also limited to Euclidean space and 
we use the two methods described above to explore the 
possible shape of the spectral function. The truncation 
scheme is detailed in section [XT] and has been shown to 
reproduce the chiral and deconfinement transition tem- 
peratures of quenched lattice QCD in Ref. [H} . In section 
IIII we briefly present our results for the quark dressing 
functions in Euclidean space below and above the critical 
temperature T c . We then introduce the spectral repre- 
sentation for the quark propagator in section IIV1 Using 
fit ansatze for the spectral function, we focus on zero 
momentum and finite quark masses in section [V] before 
moving to finite momentum in the chirally restored phase 
in section PVTl The obtained results are confronted with 
recent quenched lattice QCD results from Ref. [25|. Fi- 
nally we evaluate our truncation scheme using MEM in 
section IVTa before we conclude in section [Villi and out- 
line possible future directions. 



II. QUARK DYSON-SCHWINGER EQUATION 
AND EMPLOYED TRUNCATIONS 



A. Quark Dyson-Schwinger equation 

The renormalized quark Dyson-Schwinger equation in 
the Matsubara formalism is given by 

S^ 1 {iuj pi p) = Z 2 So 1 (iuj p ,p) - T,(iu> p ,p) , (1) 

with the inverse full quark propagator S~ 1 (iu> p ,p), the 
quark wave function rcnormalization constant Z 2 , the 
inverse bare quark propagator S , t j" 1 (iwp, p) and the to- 
be-specified self-energy Y,(iuj p ,p). With the only dif- 
ference of explicitly using imaginary arguments for the 
energy in all functions, we follow the conventions from 
Ref. |32|. For the fermionic Matsubara frequencies we 



have 



(2n p + l)nT with temperature T. For the 



Dirac structure of the inverse propagators we have 

S^ x {iu p ,p) = 174 L0 p + ry -p + Z rn m(n) , (2) 
S~ 1 (iu) p ,p) = ij4U p C{iuj p , \p\) + ij ■ pA(iuj p , \p\) 

+ B(iu) p ,\p\), (3) 

i.e. the inverse full propagator is parameterized by the 
dressing functions A(iuj p , \p\), B{iio pi \p\) and C(iuj p , \p\). 
The quark mass renormalization constant Z m and the 
renormalized mass m(/z 2 ) at renormalization scale /1, to- 
gether with Z 2 , are then determined within a MOM- 
renormalization scheme. 



In Landau gauge, the self-energy takes the form 

d 3 q 



l^S{iu) q ,q) 



x T v (iu q ,q,iuj p ,p)D ill/ {iuj p - iu q ,p- q) 



(4) 



with gluon propagator D^, quark-gluon vertex T„, 
gauge coupling a(fj,) and ghost renormalization constant 
Z3 . Assuming we know , T„ and the renormalization 
procedure, the DSE then translates into coupled non- 
linear integral equations for the desired dressing func- 
tions. 



B. Truncation setup 

The gluon propagator being a major ingredient in the 
quark DSE, can in principle be obtained in quenched ap- 
proximation from corresponding DSE's [33l436| . How- 
ever, this turns out to be a formidable task and the re- 
sults are still on a qualitative level only. We therefore 
use results from lattice gauge theory as input. These 
have been made available recently on a fine temperature 
grid over a temperature range from to 2.2 T c [32j |. 

In Landau gauge the gluon propagator is given 
through transverse and longitudinal dressing functions 
ZxiiiOk, |fe|) and Z-^iiGj^, |fe|), respectively, via 



Z T (iu k , |fc|) 
D^iiuk, fc) = -J . ^ 



p Uk)+ ^m P L Ak) 



(5) 



Hereby we introduced the projectors 



Pj„(fc) = (l-<W(i-M (<W-^r) 



(G) 



the bosonic Matsubara frequency il>k = inn^T and the 
short-hand notation fc M = (<Dfc, fe)^, for which k 2 = ujf. + 
fc 2 . As in the case of quark propagator and quark-gluon 
vertex we suppress color indices due to the remaining 
global color symmetry. 

The gluon dressing functions, Zt and Zl, have been 
calculated in quenched lattice gauge theory for gauge 
group SU(2) and SU(3) [Hj]. I R the quark Dyson- 
Schwinger equation we also need to evaluate the gluon 
propagator for momenta not identically to the ones of 
the lattice calculation. Therefore we use temperature de- 
pendent fits to these data. In Fig. [1] we show results for 
gauge group SU(3) which is the relevant gauge group in 
this work. The figure presents the lattice results denoted 
by data points together with corresponding fit functions 
represented by straight lines for a variety of tempera- 
tures. The explicit expressions for the fit functions and 
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longitudinal gluon dressing function transverse gluon dressing function 




p [GeV] p [GeV] 



FIG. 1: Longitudinal (left) and transverse (right) temperature dependent gluon dressing function for the lowest Matusbara 
frequency from lattice calculations compared with our fits. Shown are results for temperatures below the critical temperature 
(0.99 T c ) and above the critical temperature (f.lT c , 1.81 T c , 2.2 T c ). 



a more detailed discussion concerning the gluon propa- 
gator can be found in Ref. [32[ and shall not be repeated 
here for brevity. 

The details of the second ingredient for a closed quark 
DSE, namely the quark-gluon vertex, are yet to be ex- 
plored. First exploratory results on the mass and mo- 
mentum dependence of the vertex at zero temperature 
have been reported from lattice calculations and Dyson- 
Schwinger equations, see J3?i Hll and Refs. therein. How- 



ever, not much is known about its temperature depen- 
dence. In such a situation a viable strategy is to use 
phenomenological model ansatze for the vertex which 
are then justified by comparing results with other ap- 
proaches. This strategy has been successful in previous 
works and will therefore also be adopted here. 

In the following we employ the temperature dependent 
ansatz for the vertex used in Refs. (29|, [3(3, [H| . It is given 
by 



T v (iujq,q,iu) p ,p) = Z 3 



C(iLj q ,\q\) + C(iu> p ,\p\) | n x n A(iu q ,\g\) + A(iu p ,\p\) 
bivli ^ 1" I 1 ~ <W)7^ 



<h 



d 2 + k 2 A 2 + k 2 



Poafa) ln[fc 2 /A 2 + l] 

47T 



26^ 



(7) 



with the gluon momentum k^ = (uj p — 0J q ,p — q)^ and 
/?o = 11 for quenched QCD. The anomalous dimension of 
the vertex 5 accounts for the correct perturbative running 
coupling. It is given by S = —9/44 and we renormalize 
at a(fi) = 0.3. The Dirac structure with the quark vec- 
tor dressing functions A, C represents the first compo- 
nent of the Ball-Chiu vertex at finite temperature. It is 
motivated by the Slavnov- Taylor identity for the vertex 
and constitutes a first approach towards a more realistic 
non-trivial temperature dependence of the vertex. In the 
additional phenomenological fit function we have param- 
eters d\ and d 2 and the temperature independent scale 
A = 1.4 GeV. Since we use lattice results as reference 
for the gluon propagator, the quark-gluon vertex is the 



only uncertainty source in our calculation. We adjust 
the parameter d\ temperature dependent while keeping 
d 2 — 0.5 GeV 2 independent of temperature. As for d\ it 
turns out that its temperature dependence is in one-to- 
one correlation to the thermal masses mj- of the quarks as 
extracted in scctionfVl Since the same quantity has been 
determined from quenched QCD results for the Landau 
gauge quark propagator, see Tab. II of Ref. J2f|, we ad- 
just d\ such that it roughly coincides with these results. 
This procedure works only above the critical tempera- 
ture where thermal masses can be extracted. Below the 
critical temperature no such information is available and 
we therefore choose the constant value d\ = 4.6 GeV 2 . 
This choice ensures spontaneous chiral symmetry break- 



ing and agrees with the one used in Ref. [32[. We cer- 
tainly checked the d\ dependence of our results for the 
quark spectral function below the critical temperatures; 
this is discussed in section [V] Our values of d\ and the 
obtained thermal masses are summarized in Tab.|TJ Obvi- 
ously, the strength of the vertex ansatz at low momenta 
is significantly reduced above the critical temperature. 
This is what one may expect: In the Dyson- Schwingcr 
equation for the quark-gluon vertex, discussed in detail 
in Ref. [HI , we find skeleton diagrams containing dressed 
gluon propagators. As is evident from Fig. Q] there is a 
rapid decrease in the strength of the longitudinal gluon 
at or around the critical temperature T c . This rapid de- 
crease will backfeed into the vertex DSE which probably 
explains the reduction of the strength of the quark-gluon 
vertex. This behavior fits nicely to our findings for the 
parameter d\ in our ansatz for the vertex. We therefore 
believe that our vertex ansatz accurately reflects quali- 
tative and maybe even quantitative features of the fully 
dressed quark-gluon vertex. 



III. MATSUBARA PROPAGATOR BELOW 
AND ABOVE T c 

Before we focus on the spectral functions we briefly dis- 
cuss our results for the Matsubara propagator in momen- 
tum space shown in Fig. [2] The results are obtained em- 
ploying a sharp cutoff A for the computation of the self- 
energy. Thus the integration and summation extends to 
momenta and frequencies with uj 2 + q 2 < A 2 . For \n q \ < 
39 the summation over Matsubara frequencies is per- 
formed explicitly. The remaining sum is approximated 
by an integral and we checked that the results are insen- 
sitive to a change of the number of explicitly summed 
Matsubara modes. We then use a MOM-rcnormalization 
scheme with renormalization conditions C(iu!o,fi) = 1 
and B(iujQ,fi) — to(/z) and w 2 + fj, 2 = 10545 GeV 2 ; con- 
sequently all dressing functions are independent of the 
cutoff A. More details concerning the numerical proce- 
dure can be found in Ref. [30l ] . 

In Fig.[5]we show results for m(/x) = 25 MeV. Above T c 
the interaction strength as parameterized by the quark- 
gluon vertex is strongly reduced and the dressing func- 
tions A(iuio,p) and C(iojQ,p), which are finite in Landau 
gauge, do not deviate much from unity. In the infrared 
all dressing functions are constant and we observe that 
A(iuio,p) < C(iuJo,p). This is also found for the dress- 
ing functions in HTL approximation . where this in- 
equality holds for all momenta. In the simplified case 



TABLE I: Temperature dependent parameter for the quark- 
gluon vertex and extracted thermal masses tut- 
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FIG. 2: Momentum dependence of the dressing functions 
A, C, B below and above the critical temperature. Results 
are shown for the lowest Matsubara frequency and finite bare 
quark mass m(fi) = 25 MeV. 



of constant dressing functions this is perfectly consistent 
with our expectations: As can be seen from Eq. ([3]), we 
would interpret y/A/C = v < 1 as the velocity of the 
quasiparticles, \Jv 2 p 2 + B 2 /C 2 as their dispersion rela- 
tion and 1/C < 1 as the wave-function renormalization 
of the 74-componcnt 1 . In contrast to the HTL approxi- 
mation there is, however, a momentum range (8 GeV 2 to 
200 GeV 2 in Fig. O where A(iu ,p) > C(iu ,p) for all 
here considered temperatures above T c and all masses. 
To what extent this behavior is generic or truncation- 
affected needs to be investigated in future studies. 

Below T c , on the contrary, the results become less in- 
tuitive. Due to the stronger interaction and lower tem- 
perature, the chiral symmetry breaking dressing function 
B is significantly larger in the infrared with its size re- 
lated to the scale of dynamical symmetry breaking. For 
the dressing functions A(iujQ,p) and C(iu}Q,p) we find, 



1 Note that o nly the 74-component of the quark propagator obeys 
a sum rule 12311. 
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however, A(iu>o,p) > C(iu>o,p). This does not allow for 
a simple quasiparticle interpretation as outlined above. 
Also the contribution from the self-energy is of the same 
order as the free propagator, in particular the dressing 
function A(iu)Q,p) is strongly enhanced in the infrared. 
This also hints towards more complicated dynamics than 
a gas of quasiparticles, which is of course expected in the 
confined phase of QCD. 

From these results we therefore conclude that a simple 
interpretation in terms of the excitation spectrum seems 
not appropriate below T c . Even for the analysis above 
T c , as we discuss in the following, the quasiparticles can 
receive a thermal mass in the chirally restored phase. 
This is already known from the HTL expansion, where 
the quasiparticles in the high temperature and small cou- 
pling regime have a mass of order gT. 



IV. QUARK SPECTRAL FUNCTIONS AND 
REPRESENTATION 



Causality and resulting analyticity imply that a corre- 
lator can be represented through a spectral function. For 
the quark propagator this representation is given by 



duj' p(uj',p) 



2ir ioj n — lo' 



(8) 



where the Dirac structure of the spectral function is pa- 
rameterized as 

p(w,p) = 27r(p 4 (o;, |j»[)-y 4 + Pv(w, \p\) (if • P) /\p\ 

-A(w,|p|)). (9) 
Our conventions are chosen such that the scalar dress- 
ing functions themselves agree with those introduced in 
Ref. [2|| using Minkowski space conventions 2 . Assuming 
a positive definite Fock space 3 , the dressing functions 
furthermore obey 



P4(w, \p\) > vVv(", |p|) 2 + Ps (cj, \p\) 2 > (10) 
and the sum rules 



1 = Z 2 J duj pi(ui,\p\) . 

dujp v (uj, \p\) 



dujp s (uj, \p\) , 



(11) 
(12) 

(13) 



where here Z 2 is the wave function renormalization con- 
stant and not the plasmino residue which will be intro- 
duced later. 

In the following we limit ourself to vanishing momenta 
\p\ and the chirally restored phase, respectively. For this 
purpose it is instructive to introduce the projectors 



P±(P) 



(1 =R747'P/|P|) 
(IT 74) , 



(14) 



where the signs are again chosen in order to agree with 
common Minkowski space conventions. P±(p) can be 
interpreted as energy projectors for massless modes. 
For vanishing momentum it is then convenient to write 

p(u,0) =p? ML +7 4 + p M M£-74, (15) 
S(iu n , 0) = Si 1 (iuj n )L +li + SM(iu> n )L-~, 4 , (16) 

i.e. p±{uj) = 2ir(p 4 (uj, 0) ± p B (u,0)) and p y (u,0) = 0. 
From Eqs. (jTU)) and (TT3]) we infer that the functions 
p±(u>) are positive semi-definite and normalized to 2ttZ2- 
Furthermore the spectral representation yields the scalar 
relation 



duo 1 pf{u') 



2tt iui„ 



(17) 



On the other hand for a chirally symmetric phase with 
B(iw, \p\) = and p s (uj, \p\) = 0, we introduce 



p(LJ,p) = J2pe("Ap\)Pe(p)74, (18) 

e=± 

S(iw n ,p) = J2 S' e p (a;,|p|)Pe(p)74, (19) 



e=± 



i.e. p±(u, \p\) = 27r(p4(a;, \p\)±p v (ui, \p\)). As before we 
see that p±(ui, \p\) is positive semi-definite and normal- 
ized to 27rZ2 . The spectral representation for the dressing 
functions again take the form 



S£(iWn,\p\) 



du>' p±(u', \p 
2ir 



lUJ r , 



(20) 



In the following we aim to invert the linear relation be- 
tween Matsubara propagators S± (iuj n ), S±(icj n ,p) de- 
termined by solving the truncated DSE and the respec- 
tive spectral functions p± (u) and p±(w, \p\). Strictly 
speaking, for a finite set of Matsubara frequencies, this 
problem is however ill-posed. 



QUARK SPECTRAL FUNCTIONS AT ZERO 
MOMENTUM 



2 Tm ' n common Minkowski space conventions are related to those 
used in this work via 7P, = —74 and = — 17^. 

3 For completeness we note that this is not guaranteed for a gauge- 
fixed Yang-Mills theory. 



In this section we focus our attention to the quark 
propagator at zero momentum. We analyze the corre- 
lator at temperatures below and above the critical tem- 
perature and for various bare quark masses. Since the 
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general problem of extracting the spectral function from 
the spectral representation is ill-posed, we will limit our- 
selves to parameterized ansatze and determine the best 
fitting function in the given subspace. Above the criti- 
cal temperature we study the quark mass dependence of 
the so obtained insight on quasiparticle excitations. This 
strategy was explored in Refs.[24|, HBI in the framework 
of lattice QCD. 

We will mainly consider the following two-pole ansatz 
for the spectral function 



p^{uj) = 2it[Z 1 8{uj T E x ) + Z 2 5{uj ± E 2 ) 



(21) 



with fitting parameters Z\, Z 2 and E±, E%. This will be 
related to our quark propagator via Eq. (|17[) . for which 
we identify 



C DSE, A//- \ 



iuj n C(iu n ,0) ± B(icj„, 0) 
'u 2 l C 2 (iu; n ,0)+B 2 (iio n ,0) 



(22) 



Our choice is suggested by HTL results at high tempera- 
tures and small coupling. In contrast to a non-interacting 
fermion, whose spectral function consists of a single 6- 
function, it is known from HTL that an additional col- 
lective excitation develops: the plasmino. The fitting pa- 
rameters Ei t 2 denote the quasiparticle energies and Zi t2 
the corresponding residues. Certainly the full spectral 
function will be more complicated. The working assump- 
tion here is, that if quasiparticle excitations with small 
decay widths exist, then their peaks will be the domi- 
nant contribution of the spectral function and therefore 
such a simple ansatz may reveal characteristics of the 
quasiparticles and their dispersion relations. However, 
we emphasize again that in HTL calculations the ther- 
mal masses E\i<z are of order gT, whereas the width of 
the corresponding peaks in the spectral functions are of 
order g 2 T. In a strong coupling regime this becomes of 
the same order, if not larger. 

Likewise as done in Ref. [25| we also investigated other 
ansatze for the spectral function. We considered a single 
pole ansatz 



jO± (w) =2ttZ 1 5(ujtE 1 ), 



(23) 



and ansatze allowing for one- and two-particle excitations 
with Gaussian widths 



i.e. we implicitly assume uncorrelated data with equal 
total errors. For comparing the different fit forms we 
with N^j = 39. In principle, one may also 



evaluate l\ 



define £±/doi where dof denotes N u minus the number 
of fit parameters however this does not affect our results. 

A remark on the numerical results seems appropriate 
here. Within a given truncation scheme the error in 
Dyson-Schwinger calculations is essentially determined 
by the error of the numerical integration. It can be 
estimated by varying the numerical parameters of the 
integration and can be made significantly smaller than 
in lattice calculations. At least for large temperatures 
where the Matsubara frequencies are largely separated, 
the assumption of uncorrelated errors seems reasonable. 
In principle, the self-energy is calculated in DSE calcu- 
lations and the total error appears for the inverse prop- 
agator, which grows linearly in frequencies u) n . By using 
Eq. (|26|) we therefore enhance the importance of small 
ui n . We expect the relevant information for spectral func- 
tions to be encoded here as the propagator approaches 
the pcrturbative result at large frequencies very fast. 

For the data analyzed we find t\ more than one order 
of magnitude smaller for fits based on Eq. (f2"Tj) than for 
the single pole Ansatz Eq. (|23l). This is in accordance 
with the findings in Ref. [25(. Comparing obtained 
from the single peak ansatz Eq. (|24|) with its value ob- 
tained for ansatz (|21j) we find similar results for light 
quarks {m/T < 0.2) but roughly one order of magnitude 
difference in favor of (|2~lj) for heavy quarks {m/T > 0.2). 

Examining ansatz (|25j) we find parameters with im- 
proved t\ compared to the two pole ansatz. However 
depending on the initial parameter values different local 
minima may be found. Typical solutions correspond to 
a normal mode with a width which tends to zero and a 
plasmino mode with broad width. Furthermore we find 
that the possible excitations do not acquire a thermal 
mass in the chiral limit. In Ref. [25| it was found that 
ansatz (|25|) may improve x 2 /dof for uncorrelated fits but 
reduces to two delta functions for correlated fits. By way 
of comparison and taking into account the results of [25[ 
for correlated data we focus the discussion to the two- 
pole ansatz in the following. 



pfH=2^F- F i exp. 

J- 1 



>T£i) 

r? 



Z x -(^t£i) : 
-f^exp- 
t i 



Z 2 

-f^exp- 

t 2 



-{oo±E 2 f 

1 2 



(24) 



(25) 



Here Ti and T 2 are additional fitting parameters. 

In order to evaluate the quality of the fit we minimize 



'•I 



E C DSE, M I ■ s 
\S± [lUn) 



(26) 



A. Results above and below the critical 
temperature 

Before analyzing the quark mass dependence of the 
fit parameters we investigate the propagator below and 
above the critical temperature. Here, T c is given by the 
deconfinement transition as extracted from the lattice 
simulations of our input gluon propagator and from the 
dressed Polyakov loop extracted from the quark propaga- 
tor, see [32l | for details. The chiral transition temperature 
coincides with T c within a few MeV, depending slightly 
on the bare quark mass. 

In Fig. [3] we show the Schwinger function defined by 



the Fourier transform of (iu n ): 



Sl I (T) = -Tj2e~ l ^ T S I l I (iu n ) 



(27) 



The numerical data are denoted by data points. Above 
T c the solid lines are obtained by using the spectral func- 
tion obtained from fitting with ansatz (|21[) . Note that 
the points close to tT = and tT = 1 are subject to 
boundary effects due to the finite number of Matsubara 
modes used when Fourier transforming. 

At zero temperature the Schwinger function has been 
used to search for positivity violations in the quark 
propagator. From the Osterwalder-Schrader axioms of 
Euclidean quantum field theory one knows that the 
Schwinger function needs to be strictly positive for phys- 
ical particles. Conversely, positivity violations in the 
Schwinger function arc a sufficient criterion for the ab- 
sence of the corresponding particle from the physical 
spectrum of a theory. Using the spectral representation 
(jTTJ) in Eq. (|27p it is also easy to see, that a positive 
spectral function always gives a positive Schwinger func- 
tion. For the quark propagator at zero temperature these 
properties have been investigated in a number of works, 
see [4(| and Refs. therein. At finite temperature it has 
been suggested from calculations in a simple model, that 
positivity violations only occur below the critical tem- 
perature, whereas positivity is restored above T c [41 [. 



Although this is not exactly what we see in our more 
elaborate calculation, there are clear signals for a qual- 
itative change in the Schwinger function at T c . This is 
in line with quenched lattice QCD results that directly 
determine S^(t) HI). 

Above the critical temperature S+(t) is positive and 
found to be convex on a log-scale. Furthermore chiral 
symmetry restoration through B(iw n ,0) = translates 
into Sf(ico n ) = -S^(-i(J n ) and S^(t) = Sjjf (1/T-t). 
We see this symmetry emerge when decreasing the cur- 
rent quark mass m(/z). Related to dynamical symme- 
try breaking this symmetry is absent below T c , even in 
the chiral limit. More striking, however, is the concave 
shape around tT « 0.5 on a log-scale. This cannot be 
reproduced by any ansatz including real poles only. This 
finding again agrees with quenched lattice calculations in 
Ref. [2j|. 

In addition we find positivity violations at T = 0.99 T c 
for a range in current quark masses, albeit this effect is 
not seen in the chiral limit and also for the smaller tem- 
perature. Nevertheless, the concaveness of the Schwinger 
function below T c is generic in the sense that it is in- 
dependent of variations in d\. It can also not be at- 
tributed to the dynamical generation of quark masses: we 
find concave curvature below the deconfinement transi- 
tion temperature T c also in the chiral limit and for values 
of d\ that do not allow for spontaneous dynamical mass 
generation. It is therefore apparent, that the curvature 
of the Schwinger function is related to quark confinement 
[4l| and may serve as another source for extracting the 
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FIG. 3: Schwinger function (t) for different temperatures 
and current quark masses. For T — 0.99T C and m(^i) = 
0.2, 0.4 GeV the Schwinger function changes sign in an in- 
terval around tT m 0.8. 



deconfinement transition temperature T c from the quark 
propagator besides the dressed Polyakov loop [29l. |3U |42| . 



B. Quark mass dependence of fit parameters 

In the upper part of Fig. H] the dependence of the poles 
Ei, E-2 on the current quark mass for three different 
temperatures is shown. The lower part shows the rel- 
ative strength of the plasmino pole, defined as the E2 
branch in the upper plot, which is given by the ratio 
Z2/ {Zi + Z2 ) • The disappearance of the scalar part of the 
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T/T = 2.2 



T/T c = 1.5 
T/T = 1.25 

c 



0.25 0.5 0.75 1 1.25 
m(n)/T 



FIG. 4: Dependence of fit parameters in two-pole ansatz at 
vanishing momentum on current quark mass m(/i) for three 
temperatures. For the two branches in the upper plot we 
define E 2 >E 1 . 



propagator in the chiral limit, i.e. B(iuj n ,0) — 0, mani- 
fests itself in the spectral function being an even function, 
p(—uj, 0) = p(cj, 0), and therefore E\ = E%, Z\ = Z%. We 
use this as the definition of the thermal mass mj = E\ 
and, as outlined in the context of Tab. H use this value 
to determine the vertex parameter d\ . It is worth noting 
that the qualitative behavior in Fig. @] does not change 
when increasing d\ by one order of magnitude. 



From the mass dependence of Z 2 /(Zi + Z 2 ) in Fig. @] 
it is discernible that the plasmino contribution decreases 
with increasing quark mass. This is in accordance with 
the expectation that for heavy quarks the spectral func- 
tion reduces to one of free quarks. Only small deviations 
from sum rule (jllj) are found with a mismatch of less 
than 4%. Furthermore we find the minimum of Ei at 
non- vanishing quark mass whereas the plasmino pole in- 
creases monotonically with quark mass. The slope of Ei 
and E2 for fixed m(fi)/T is temperature dependent and 
decreases with increasing temperature. Compared to the 
lattice data in Ref. [2f| we find qualitative and quanti- 
tative similar results. The minima in E\ approximately 
coincide with the lattice results whereas the slopes in 
Z 2 /{Zi + Z 2 ) and Ex 2 are slightly different. However 
this may be traced back to different definitions of the 
bare quark mass. 



VI. MOMENTUM DEPENDENCE OF THE 
SPECTRAL FUNCTION IN THE CHIRAL LIMIT 

In this section we investigate the quark propagator in 
the chiral limit for varying momenta. For a chirally sym- 
metric quark propagator it is convenient to use the de- 
composition on energy projectors as shown in Eqs. (|18[) 
and (|19p . In terms of the quark dressing functions we 
then have 



S£(iuj n ,\p\) 



iu n C(iu n , \p\) ± \p\A(iuj n , |p|) 
4C 2 ( lWn ,|p|) + |pPA2^„>[ 



(28) 



Along the same line as in the previous section we use the 
two-pole ansatz 

p£(w, \p\) = 2k[Z 1 5{uj TE X ) + Z 2 8{u ± E 2 )] (29) 

and fit the parameters Z\, Z 2 and E\, E 2 to the data. 
These are then momentum dependent and we obtain the 
dispersion relation of quasiparticle and plasmino together 
with there spectral strength. For all momenta shown here 
i 2 ^ is of the same order of magnitude than the results for 
zero momentum. 

In Fig. [5] we show the results for the momentum de- 
pendence of the parameters E\, E 2 and Z 2 /(Zi + Z 2 ) for 
temperatures T/T c = 1.25, 1.5 and 2.2. For increasing 
momenta the ratio Z 2 j[Z\ + Z 2 ) decreases and tends to 
go to zero for high momenta. Likewise E\ is approaching 
the light cone. Therefore the spectral function at high 
momenta reduces to the one of a single free particle. In 
the plasmino branch a minimum at non- vanishing mo- 
mentum is clearly recognizable. This feature is known 
from perturbation theory and may be anticipated from 
the results in Ref. [25|. It is conspicuous and not ob- 
tained in one loop perturbation theory that the plasmino 
branch goes over into the space-like region. However it 
is known that due to Landau damping the spectral func- 
tion at high temperatures contains a contribution in the 
space-like region which is not included in our fit ansatz 
Eq. (|29p . It might be that the plasmino pole by shifting 
into the space-like region mimics this missing contribu- 
tion. This possibility is also considered in the following. 
An analytic continuation of the DSE to include smaller 
Euclidean energies might help to distinguish these differ- 
ent scenarios. Also note that with increasing momenta, 
the value of E 2 becomes less constrained by the fitting 
procedure, since the spectral strength Z 2 and therefore 
the contribution of the plasmino branch to the data de- 
creases. 

In Fig. [6] we focus on the region near zero momentum. 
In addition to the fitted values shown as straight lines, 
we also present the behavior for the HTL result [39( given 
by 
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FIG. 5: Momentum dependence of the fit parameters for the 
two pole ansatz (|29l) at three different temperatures. 



FIG. 6: Close up of Fig.[S]in momentum range p/T E [0, 0.2]. 
For comparison we give the small momentum behavior of the 
HTL approximation in Eqs. (|30l) and (|3I|) as dashed lines. 



for |p| <C rriT- We observe that our results agree well with 
the HTL result for sufficiently small momentum, with the 
better agreement for the plasmino branch E 2 , and then 
depart from it. Since the expansion is supposed to apply 
for small |p|/T, this is expected. Comparing the two 
temperatures it is recognizable that the results at higher 
temperature agree better with the HTL result, although 
we are certainly not in a regime where we expect the 
HTL approximation to work. 

Since the result for the plasmino branch of the two-pole 
ansatz indicate a possible space-like continuum contribu- 
tion we extend the two-pole ansatz (|29|) by the HTL con- 
tinuum contribution. The fitting function then matches 
exactly the spectral function in HTL approximation 

P±(w, \p\) = 2tt[Zi*(w T Ex) + Z 2 S(lo ± E 2 )] 

+ ^- m 2 T {1 T - x 2 ) 

\P\ 



(\ P \(lTx)±^[(lTx) 
2 



x In 



X + 1 


± 2 




x- 1 


) 







(32) 



where x = u/\p\- The thermal masses are taken from 
Tab. U and Ex 2 and Z x 2 are the fit parameters. The 
results are shown in Fig. [7J We find i 2 + to be marginally 
smaller compared to the two-pole ansatz and we conclude 
from the results presented in Fig. [5l that the plasmino 
branch going over into the space-like region might well 
be an artifact of a missing space-like spectral weight in 
ansatz (f2Tj) . The upper plot of Fig. [JJ shows that the 
plasmino tends to the light cone for momenta p/T ~ 0.8 
instead of crossing it. Also the minimum of the plasmino 
branch has shifted to lower momenta. At high momenta 



the plasmino seems to deviate from the light cone again. 
However we note again that due to the small spectral 
strength Z 2 at these momenta the value of E 2 is only 
poorly constrained. On the other hand the momentum 
dependence of Ex is practically unaltered. 

The lower plot of Fig. [JJ shows the relative strength of 
the plasmino Z 2 /Z toi where 



+ |p| 



Z t , 



Z x +Z 2 



dujp^(u, \p\ 



(33) 



-IpI 



Compared to the two-pole ansatz the plasmino spectral 
weight decreases considerably more rapid with increasing 
momentum. We also find that Z to t agrees within 3% with 
the expected value obtained from sum rule (fTTj) . Since 
the spacelike part is sublcading for small momenta the 
small momentum behavior of Ex %2 and Z 2 j[Zx + Z 2 ) is 
practically unaltered by the continuum contribution. 



VII. QUARK SPECTRAL FUNCTIONS VIA 
MEM 



In order to shed more light on the possible shape of the 
spectral functions and also on the truncation for gluon 
propagator and quark-gluon vertex 4 , we employ MEM. 
For simplicity our discussion will be based on historical 



4 Note that their parameterization used here is only constrained 
in Euclidean space, but their analytic continuation is in principle 
relevant for the shape of the quark spectral function. 
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FIG. 7: Momentum dependence of the fit parameters for the 
extended ansatz (|32p for two different temperatures. 



MEM (see e.g. 
functional 



Rcfs. [lj| Hll), where we maximize the 



Q[p]= a S[p]-- X 2 [p] 



(34) 



in the spectral function p. Its constituents are the en- 
tropy S[p) and the % 2 [p]-distribution. In the case of a 
scalar channel as in Eq. (|17|) and Eq. ([20|) it is common 
to use the Shannon-Jaynes entropy 



S[pj 



P± ( w ) — m(uj) — p± (w) log 



m{ui) 



(35) 



for the entropy functional with prior estimate m(w). Fea- 
tures of this choice are that there is only one minima 
in Q[p] for each choice of m(u>) and that p(ui) has the 
same sign as m(w). For the x 2 [/^-distribution between 
the propagator given through p(oj) and the data from 
DSE calculations we will assume 



x 2 [p±] = \t 



(36) 



As discussed before this corresponds to assuming uncor- 
rclatcd data with equal errors and we will need to specify 
the variance a 2 . For the prior estimate, which essentially 
determines the support of the spectral function and its 
semi-positivity, we assume 



m(w) = m #(A 2 - w 2 ) 



(37) 



The magnitude tuq and the cutoff A will be varied in the 
following to analyze their influence on the shape of the 
extracted spectral function. The method as presented 
here follows the philosophy of regularization [19] : Since 
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the inversion of the spectral representation is an ill-posed 
problem, we want to impose the positivity condition on 
the spectral function in a way that leads to a well-defined 
problem. 

In historical MEM the Lagrange multiplier a is then 
chosen such that x 2 [p±] — N u . Within our setup we 
observe, however, that the extracted spectral functions 
p±(u>) actually only depends on the combination a = 
cut 2 and that in this case a 2 = £±/N u . We will exploit 
this feature by varying a in the following, in particular 
because we do not want to give an estimate for a 2 at this 
point. 

In order to get an intuition for the way MEM works, we 
consider p — and produce 'mock data' by determining 
S+(tj„) for the spectral function 

p + (w)=27rQtf(M-a;) + ^(M + w)) (38) 

analytically. This corresponds to the two-pole 
ansatz (p?Tj) and also with the HTL result at p = 0. 
Our results are summarized in Fig. [5J By the example 



11 



10"' 



10 



> 
o 



to 



10" 



10" 



! | Hill 

— 


| lllllll 


| lllllll 


linn 


I linn 


| lllllll 


1 












Mr 
































- 























m = 10" 2 GeV" 
m = 10" 2 GeV" 
m = 10" 2 GeV" 


, A= 2 GeV 
, A=5 GeV 
, A= 10 GeV 























m = 10~ 2 GeV~ 


, A= 20 GeV 












m = 1 GeV 


, A= 2 GeV 




- 1 ■ 


1 


1 ■ 




m = 10~" GeV 


', A=2GeV 




■ 


1 ■ 


1 ■ 





10" 



10" 



10_ 10- 
ot [GeV 2 ] 



10" 



10 



10" 



FIG. 9: l\ between numerical DSE results and S+(ui n ) for 
MEM extracted spectral function as a function of a for various 
mo and A as well as T — 1.5 T c and = 20. 



of mo = 10~ 2 GcV _1 and A = 2 GeV, we observe that 
the correct spectral function is reproduced when making 
a sufficiently small. This is very generic, in particular we 
checked that the extracted spectral function for a given 
a G [l(r 15 GcV" 2 ,lCT 5 GcV" 2 ] is indistinguishable on 
the plot when varying mo G [1 GeV -1 , 10 _4 GeV _1 ] and 
A G [2 GeV, 20 GeV]. For the corresponding t\, which 
should be equal to N u a 2 by choosing an appropriate 
a, we find on the one hand side, that the dependence 
on the prior estimate is most extreme at larger a, al- 
beit in general the sensitivity can be considered as mild. 
More importantly, we find that l\ for different prior es- 
timates converges quickly upon decreasing a and that it 
converges towards zero. This can be qualitatively under- 
stood from the fact that we know that there is a positive 
definite spectral function which reproduces the data and 
that there is a competition between a5[/9+] and ^x 2 [/>+] 
when maximizing the functional Q[p+]. This leads to 

> const > 0. Conse- 



\x 2 [p+] and S[p^ 



aS[p+] 

qucntly, v\ should go to zero when taking a — > 0. 

Before discussing errors on our computed DSE 'data', 
we want to present the results when applying MEM to 
our results in a similar fashion as discussed above. For 
this purpose we present (. 2 + as a function of a at p = 
for various prior estimates in Fig. [9] and illustrate some 
extracted spectral functions for too = lCP 2 GeV -1 as well 
as A = 10 GeV and A = 20 GeV in Fig. [TUJ 

The two-pole fit ansatz (|2Tj) for this specific case 
yields i 2 ^ ~ 10 -5 GcV~ 2 . As anticipated we obtain bet- 
ter fitting spectral functions for small enough a, here 
a < 10 _5 GeV -2 for the considered prior estimates too G 
[lGeV" 1 ,10- 4 GeV _1 ] and A G [2 GeV, 20 GeV]. As be- 
fore we also find that the dependence on mo is very weak 
and we can therefore focus on the spectral function's sup- 
port given through A only. It is, however, striking that 
the i\ has a strong dependence on A. Of course we 
should be able to obtain better fitting spectral function 
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FIG. 10: Upper panel: Spectral functions extracted from 
numerical DSE results at various values of a for mo = 
10~ 2 GeV~\ A = 10 GeV and N u = 20. Lower panel: Same 
us upper plot but for A = 20 GeV. 



when increasing the support, the results suggest however 
that a few percent of the spectral strength come from en- 
ergies |w| > 2GeV near A. This is unexpected. As shown 
in Fig. [TU] this does not affect the shape of the extracted 
spectral function in the domain ui G [—2 GeV, 2 GeV] for 
a ps 10 _8 GeV -2 noticeably. In case of the mock data we 
could already observe a two-peak structure at this point. 

When decreasing a beyond 10 _8 GeV -2 , we find a 
strong dependence on A and even more puzzling a fast 
saturation of l 2 + . The fact that the spectral functions 
tend towards a sum of peaked functions should be con- 
sidered a generic feature of method and not be over- 
interpreted. The agreement at A = 20GeV with the two- 
pole ansatz is probably a coincidence, in particular since 
the true spectral function is expected to have a finite 
width. 

The reason for our conservatism is the value of l\ 
that can be obtained within MEM. Even for an al- 
lowed support of [— 20GeV, 20GeV], we cannot obtain 



£ 2 + < 8 10~ 10 GcV On the other hand a variation 
of numerical parameters suggests that the numerical er- 
ror on our 'data' is in fact below this value. Therefore, 
strictly speaking, historical MEM requiring \ 2 = N u can- 
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not be applied. 

A plausible explanation for this finding would be a 
spectral function that, at least for the applied trunca- 
tion scheme, is not positive semi-definite. To elaborate 
this a bit further, it is known e.g. from the calculation 
in Ref. [l8|, that the determination of the spectral func- 
tions via the quark's DSE requires knowledge of gluon- 
propagator and quark-gluon vertex in the time-like do- 
main and in particular of the gluon's spectral function. 
For the parameterization of the gluon used in this work, 
it can be checked that its spectral function is not positive 
semi-definite. 

We therefore consider the analysis via MEM on the one 
hand side as an advise for caution when looking at the 
fits discussed in the earlier sections. On the other hand it 
motivates future work for the direct determination of the 
spectral function and a better analysis of the truncation 
scheme. 



VIII. CONCLUSIONS AND OUTLOOK 

In the present work we have analyzed the quark prop- 
agator at finite temperatures obtained within a recently 
employed truncation of the quark's Dyson- Schwinger 
equation. The investigation was in part motivated by 
lattice QCD calculations of the quark propagator in Lan- 
dau gauge using the quenched approximation, which - in 
addition - also attempted to extract the spectral function 
of the propagator from the obtained data. 

For the Matsubara propagator, i.e. in Euclidean space, 
we found qualitative agreement with the lattice QCD re- 
sults below and above the deconfinement phase transi- 
tion. General analytic properties of the quark propagator 
are different below and above the deconfinement phase 
transition temperature and as our results indicate, de- 
termined by the non-perturbativc behavior of the gluon 
propagator. Consequently, we also conclude that our nu- 
merical data below the deconfinement phase transition do 
not allow for an interpretation in terms of a quasi-particle 
picture. Above the deconfinement phase transition, on 
the other hand, we could fit the data by a physically 
motivated spectral function consisting of a quasiparticle 
branch and a plasmino branch with vanishing width, re- 
spectively. The obtained results for the current quark 
mass dependence of the thermal masses and for the dis- 



persion relations of quasiparticle and plasmino are then 
consistent with those from quenched lattice QCD. Fur- 
thermore we included a possible spacelike continuum con- 
tribution in the model spectral function which was not 
considered in the lattice calculations. We find this con- 
tribution to shift the dispersion relation of the plasmino 
into the timelike region. This indicates the importance of 
additional continuum contributions in the spectral func- 
tion. 

In addition we also analyzed our results by means of 
the so-called historical maximum entropy method. Al- 
though the method was shown to work for 'mock' data 
once the input was assumed to be accurate enough, the 
application to the numerical DSE results turned out to 
be more subtle. We interpreted the results as a hint that 
the spectral function, at least in the employed truncation, 
might not be positive semi-definite. 

In summary, we conclude from this that we can achieve 
nice agreement with lattice QCD data in Euclidean space, 
but extracted spectral functions using fit forms should 
be taken with a grain of salt. Within the framework of 
Dyson- Schwinger equations this motivates future work of 
directly evaluating the spectral function. Although tech- 
nically more evolved, calculations within models or using 
simpler truncations indicate that this is doable. Also the 
sensitivity on the gluon propagator and the quark-gluon 
vertex should be investigated, in particular since their 
analytic continuation to Minkowski space and therefore 
their spectral function is relevant for the outlined calcula- 
tion. Furthermore it might be illustrative to analyze the 
lattice results of the gluon propagator along similar lines 
as has been done for the results of the quark propagator 
in quenched approximation. 
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